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Abstract. 

We consider dynamical systems depending on one or more real parameters, and assuming that, for some "critical" value of the 
parameters, the eigenvalues of the linear part are resonant, we discuss the existence — under suitable hypotheses - of a general 
class of bifurcating solutions in correspondence to this resonance. These bifurcating solutions include, as particular cases, the 
usual stationary and Hopf bifurcations. The main idea is to transform the given dynamical system into normal form (in the 
sense of Poincare-Dulac), and to impose that the normalizing transformation is convergent, using the convergence conditions 
in the form given by A. Bruno. Some specially interesting situations, including the cases of multiple-periodic solutions, and of 
degenerate eigenvalues in the presence of symmetry, are also discussed with some detail. 
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1. Introduction 



In this paper we consider dynamical systems of the form 

u = f{u,X) = A{X)u + F{u,X) u(t)eR" A € (1) 

depending on one or more real parameters A, and, assuming that, for some "critical" value A = Aq of the 
parameters, the matrix A(Ao) admits resonant eigenvalues, we want to discuss the existence - under suitable 
hypotheses - of a general class of "bifurcating solutions" u = u\ (t) in correspondence to this resonance. These 
bifurcating solutions include, as particular cases, the usual stationary bifurcation and the Hopf bifurcation 
as well (see, e.g. [1]). The main idea is that to transform the given dynamical system into normal form (in 
the sense of Poincare-Dulac [2-9]), and try to impose that the normalizing transformation is convergent. The 
imposition of convergence is essentially based on the application of the conditions given by Bruno [4-5], and 
leads to some prescriptions which may be fulfilled thanks to the presence of the parameters A; in this way, 
the appearance of these "resonant bifurcating" solutions can be automatically deduced. 

Among these solutions, some attention is paid to discuss some situations with special interest (also from 
the physical point of view), including the cases of multiple- frequency periodic solutions, and of degenerate 
eigenvalues in the presence of symmetry, giving examples for each situation. 



2. Basic assumptions and preliminaries 

We need some preliminary notions and results. 

Let u = u{t) e R", A e W and / = f{u, A) be a vector-valued analytic function in a neighbourhood 

of some uq and Aq (it is not restrictive to choose uq = and Aq = 0), such that /(A, 0) = for each A in a 
neighbourhood of Aq. Let us denote the linear part of / by A{X)u where 

^(A)=V„/(A,0) (2) 

and assume that, for some "critical" value of the parameters (we can assume that this value is just Ao = 0) 
the matrix 

Ao = A{0) 

is semisimple, i.e. diagonalizable. Let us notice that in the case is not semisimplc, then it could be 
uniquely decomposed into a sum of a semisimple and nilpotent part: Aq = A^^^ + A'-"-', and the foregoing 
discussion could be equally well applied to the semisimple part A^^^ , considering in particular normal forms 
with respect only to yl^") . The introduction of normal forms with respect to a non-semisimple matrix requires 
a more difficult procedure and will not be considered here (cf. [10-11]). 

Up to a linear change of coordinates (possibly after complexification of the space), we will then assume 
for convenience that the matrix Aq is diagonal, with eigenvalues cti , . . . , (t„ . The first important assumption 
is that, for the value Aq = 0, the eigenvalues exhibit a resonance, i.e. there are some non-negative integers 
rrii such that 

n n 

niiai = aj ; ^ mj > 2 (3) 

i=l i=l 

for some index j G [1, . . . , n]. 
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Together with the given dynamical system (DS) (1), we need also to consider its normal form (NF) (in 
the sense of Poincare-Dulac [2-9]), in a neighbourhood of uq = 0, Aq = 0. As well known, the idea is that of 
performing a near-identity coordinate transformation 

u ^ V — u + (p{u) (4) 

in such a way that in the new coordinates v the given DS takes its "simplest" form. To define this, consider 
the linear operator A (the "homological operator" [2-9]) defined on the space of vector-valued functions h{v) 
by 

A{h) = Aqv ■ Vh - Aoh (5) 

Writing the NF in the form 

v = g{v,X)= Anv + G{v,X) (6) 
the nonlinear terms G{v, A) are then defined by the property [2-9,11-14] 

G{v, A) e Kcr{A) (7) 

Actually, to be more precise, one should also consider, together with (1), the p equations 

A = 



(as in the usual suspension procedure), and extend the homological operator A adding to the matrix Aq the 
last p columns and p rows equal to zero, and similarly extend G{v, A) as a (n + p)— dimensional vector-valued 
function with the last p components equal to 0. It is important indeed to notice that it is essential here to 
consider the A as independent variables; in this way, in particular, one has that the bilinear terms in v and 
A are included in G{v, A); sec also Remark 1 below. 

Let us now briefly recall the following important results. The proof of these can be found e.g. in 
[8,11-14]. 

Lemma 1. Given the matrix Aq, the most general NF is given by (6) with 

G{v,X) = Y,Pi{\Pi'"))B,v where B, e C(ylo), p = p{v) G Iao (8) 

i 

where C{Ao) is the set of the matrices commuting with Aq, and is the set of the constants of motion of 
the linear system 

V = Aqv (9) 

The sum in (8) is extended to a set of independent matrices Bi, the constants of motion p(v) can be chosen 
in form of monomials (possibly fractional) and the functions (3i are series or rational functions of the Vi (see 
[13] for a detailed statement, and below, for the cases of interest for our discussion). 

It is now clear that the assumption (3) on the existence of some resonance among the eigenvalues of Aq 
ensures that there are nontrivial constants of motion of the linear problem (9), and then nontrivial terms in 
the NF (8). 

It is also well known that the coordinate transformations taking the given DS into NF is usually per- 
formed by means of recursive techniques, and that in general the sequence of these transformations is purely 
formal: indeed, only very special conditions can assure the convergence of the normalizing transformation 
(NT). 

Let us now recall the basic conditions, in the form given by Bruno, and called respectively Condition lj 
and Condition A, which ensure this convergence. The first condition is (see [4-5] for details) 
Condition u: let Uk = m.m\{q,a) — aj\ for all j = 1, . . . ,n and all n—uple of nonnegative integers Qi such 
that 1 < Qi < and {q,(j) = qiCTi ^ aj: then we require 

oo 

fe=l 
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This is a actually very weak condition, devised to control the appearance of small divisors in the series of 
NT, and generalizes the Siegel-type condition: 

\{q,a)-aj\ >e (^Qi) 

i=l 

for some e, ;/ > 0, or the much simpler condition | (g, cr) — cTj | > e > 0, for all n— uple gj such that (g, a) 7^ crj 
(see [2-5]). We explicitly assume from now on that this condition is always satisfied. 

The other one, instead, is a quite strong restriction on the form of the NF. To state this condition in 

its simplest form, let us assume that there is a straight line through the origin in the complex plane which 
contains all the eigenvalues Ui of Aq . Then the condition reads 

Condition A: there is a coordinate transformation u — > u changing f = Aqu + F to a NF g = g{v) having 
the form 

g{v) = Aqv + a{v)AQV 

where a{v) is some scalar-valued power series {with a{0) = 0). 

In the case there is no line in the complex plane which satisfies the above property, then Condition A should 
be modified [4-5] (or even weakened: for instance, if there is a straight line through the origin such that 
all the cTi lie on the same side in the complex plane with respect to this line, then the eigenvalues belong 
to a Poincare domain [2-5] and the convergence is guaranteed without any other condition); but in all our 
applications below we shall assume for the sake of definiteness that the eigenvalues are either all real or 
purely imaginary; therefore, the above formulation of Condition A is enough to cover all the cases to be 
considered. 

Remark 1. It can be useful to point out that it is essential in the present approach to use the suspension 
procedure for the parameters A, i.e. to consider A as additional variables. Indeed, let us consider, for instance, 
a simple standard Hopf-type 2— dimensional bifurcation problem with p=l parameter: 

u = Aqu-\- XI u + higher order terms where Aq = 

and / is the identity matrix. If A is kept fixed ^ 0, the eigenvalues of the linear part Aq + XI are a = X±i 
and, as a consequence, there are no (analytic nor fractional) constants of motion of the linearized problem, 
and the NF is trivially linear. Also, the NT would be convergent, as a consequence of Condition A (or - more 
simply - of the Poincare criterion [2-5]). But no bifurcation can be found in this way. Instead, considering A 
as an independent variable, the linear part of the problem is m = Aqu and now there are nontrivial constants 
of motion in the NF, i.e. the functions of r"^ = uf + u\ and A. 

Given the DS (1), it will be useful to rewrite its NF, according to Lemma 1, observing that obviously 
Aq e C{Aq) and in view of the above Condition A, in the following "splitted" form 

v = g{v,X)= Aov + a{X,p{v))AQV + Yl]l^A^^Pi'"))B3V p{v)elAo (10) 

where (hereafter) ^* is the sum extended to the matrices Bj ^ Aq Following Bruno [4-5], we can then 
say that the convergence of the NT is granted if = 0- Clearly, "convergence" stands for 

"convergence in some neighbourhood of uq = 0, Aq = 0" . 

Let us remark, incidentally, that an algorithmic implementation of the procedure for obtaining step by 
step the NF is possible (cf. e.g. [15,16]); we stress however that actually, in this paper, we shall not need 
any explicit calculation of NF's. 
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3. The 2— dimensional case. 



For the sake of simphcity we consider first of all the case of 2-dimensional DS, with some other simplifying 
assumptions; more general cases will be considered in subsequent sections. 
Theorem 1. Given a 2-dimensional DS 

u = A{X)u + F{u,X) weR^AeR (11) 

(i.e. n = 2,p= 1), assume that for Aq = the two eigenvalues ai,a2 of Aq are nonzero, of opposite sign if 
real, and satisfy a resonance relation (3), which in this case can he more conveniently written in the form of 
commensurability of the two eigenvalues; 

^ = -^ = eo (12) 

Sl S2 

where si,S2 are two positive, relatively prime, integers. Assume also that 

S2:^^ll(A) + 51^^22 (A) ^0 (13) 
dA dA A=o 

then there is a "resonant bifurcating" solution u = u\{t) of the form 

ui =c'^^ exp (si^(A)t) + ^2 exp {nid{X)t) 

m = — oo 

U2 =c(^^ exp ( - S26'(A)t) + ^ c„2 exp {n29{X)t) 

712 = —00 
n2 -S2 



(14) 



where for A — > 



si9iX)^ai ; -S20{X) ^ a2 ; c'^^^c'^^^^O 



and all terms in the two series are "higher-order terms", i.e. terms vanishing more rapidly than the two 
leading terms, and the series are convergent in some time interval. There is also an analytic constant of 
motion along this solution, which, for small X, has the form 

p(?2) = (c(i))^^(c(2))^^+h.o.t. 

Proof. Consider the (possibly non-convergent) coordinate transformation u v (4) which takes the given 
DS into NF (10). The assumption on the resonance of the eigenvalues ai and Lemma 1 imply that the DS 
in NF is expected to have the following form 

i) = g{X, v) = Aov + a{X, p)Aov + /?(A, p)Bv = (1 + a)Aov + (5Bv (15) 

where p = p{v) = vl^V2^ is a (monomial) constant of motion of the linear problem v = Aqv, B is any diagonal 
matrix independent of Aq (i.e. such that S2-B11 + S1B22 7^ 0), and A can be viewed as a constant of motion 
of the enlarged system including A = 0. We now impose the Bruno Condition A (in the form given above, 
indeed ai are either real or purely imaginary) in order to ensure the convergence of the NT: this amounts to 
imposing 

f3{X,p) = (16) 

The expression of this function is clearly not known (unless the NF itself is known), but the manifold 
defined by /? = is analytic [4], and the relevant behaviour of the first-order terms (in A) of (3{X, p) can be 
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inferred from the original DS ii = A{\)u + . . . , indeed the first-order terms (in A) of A{X)u, i.e. Aqu + 
X[dA{X)/dX) \)^^qU are not changed by the NT; therefore 



and 



a(0,0) = /3(0,0) = 



52:^^11 (A) + 51-^^22 (A) = (S2-Bii + S1B2 

dA dA A=o ^ 



9/3 
'9A 



A=0,p=0 



having also used S2CT1 + Sia2 = 0. Assumption (13) then shows that one can satisfy the condition /3(A, p) = 0, 
thanks to the implicit-function theorem, if A and p are related by a function 



With (3 = 0, a solution of (15) is then 



or, putting ^(A) = ^o(l + ce)j 



X = X{p) with A(0) = 

v{t) = exp ((1 -I- a)Aot)vo 
vi{t) =vioexp{si9{X)t) 

V2{t) =V2oexp{-S20{X)t) 



(17) 



with the constraints 



{vi{t)Y^ {v2{t)Y^ = {vioY^ {v2oY' = P&^Ao and A = A(p) 



(18) 



(19) 



and where a = a(X{p)) , ^(A) — > 9o for X,v,p — > 0. On the other hand, in the analytic manifold 
defined by /3(A, p) = the NT is convergent and this bifurcating solution corresponds to an analytic solution 
u = u\{t) of the initial problem. The original coordinates u are in fact related to the new ones v by an 
analytic transformation (the inverse of (4)) 

u = V + tp{v) 

where ip is a power series in the fi, V2, and the convergence is granted on some neighbourhood of zero, say 
It'll < Ri, \v2\ < R2- Notice now that the condition p{v) = const does not ensure, if the eigenvalues are real 
(see the examples below), that the variables are bounded for all f e R; in this case, it will be sufficient to 
choose A, together with viq, V20, small enough in order that the convergence is granted in some time interval 
T^<t< T2. 



Example 1. Let the matrix elements Aij{X) of A{X) in (11) satisfy 

Aii(A)^l ; A22{X)^-2 ■ Ai2(A) and A2i(A) ^ for A^O 

Then ai ~ l,cr2 ~ ^2, the constant of motion of the linear DS v — Aqv is p — v1v2, and, assuming 
d(2Aii + A22)/dX\x^o ^ 0, there is a resonant bifurcating solution whose leading terms have the exponential 
behaviour 

ui =c(^^ exp ((1 + a)t) + ... 
U2 =c(^) exp ( - 2(1 -I- a)t) + ... 

with a, c^^^, c^^) — > 0. To give a more explicit case, let us assume that the first terms of the DS (11) are 

A2i{X) -2 + Xj \utulj 
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with Ai2, A21 ar})itrary vanishing functions of A. In this case, the first nonhnear term is already in NF, with 
B = I = Identity (of. eq.(lO)) and one easily gets 0{X) = 1 + A/3 + . . ., a = A/3 + . . ., whereas eq. (16) 
takes the form f3 = |A + p + ... = giving then 



A = --U,U2 + ■■■ 






The next example, even if apparently similar, is actually a generalization of Theorem 1, indeed eigen- 
values of the same sign will be involved. In this case, the convergence of the NF would be guaranteed with 
no other condition, thanks to the Poincare criterion; however, we need even in this case to impose that 
Condition A is fulfilled, in order to obtain a bifurcating solution of the same form as discussed so far. 

Example 2. Let Aq = diag(l,2), and d(2Aii(A) - A22{X)) /dX\x=o ^ 0. Then now p = v1/v2 and the 
bifurcating solution is expected along a manifold of the form Xu2 oc tif + . . ., obtained as solution of the 
condition (16). The discussion is now identical to the previous cases. •O 



Corollary 1. (Hopf bifurcation) If (J\,(J2 are imaginary, ai^2 
standard "transversality condition" 

d Re a(A) 



±icoo, then condition (13) becomes the 



dA 



A=0 



7^0 



and the bifurcating solution is the usual Hopf hifurcaiion. 

Proof. The proof of Theorem 1 holds also in the case of imaginary eigenvalues. In this case, si 
vf + V2 = r"^, d{X) = i(wo + a(A(r^)) = iui{X), and a simple calculation shows indeed that 



S2 



1, P = 



d(Aii + A22) 



dA 



d Re 



A=0 



dA 



A=0 



Notice that in this case the condition p = r^ = const ensures \vi\ < r and then the convergence is true for 

all tea. • 



4. DS with dimension n > 2: reduction to a lower dimensional problem 

In this section, we will consider a general situation in which, given a n— dimensional DS with n > 2, it is 
possible to reduce the problem to a lower dimensional case. A quite simple but useful result is the following. 

Lemma 2. Consider a n— dimensional DS (n > 2), and assume that for Aq = there are r < n resonant 
eigenvalues, say ai,. . . ,ar, such that no resonance relation of the following form 

r 

^mh<Th=(Tk k = r + l,...,n (20) 

h=l 

exists. Then, the DS - once in NF - can be reduced to a r— dimensional problem putting for the remaining 
n — r variables 

Vk{t) = k = r + 1, . . . ,n 

If this r— dimensional problem admits a solution Vh{t), h = 1, . . . ,r (e.g., a bifurcating solution as in Theorem 
1 ), then the original DS admits a solution in which the n — r components Uk{t), k = r + 1, . . . ,n, are "h. o.t." 

with respect to the first r components Uh{t). 

Proof. Let us consider the NF variables Vi, i = 1, . . . ,n, and let us introduce the shorthand notation v' = 

(A! \ 

{y\,..., Vr) and v" = {Vr+i, . . . , Vn), and - correspondingly - = ( ° ^„ j with Aq = diag(cri, . . . , ar) 
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and similarly for Aq. According to Lemma 1, the matrices Bi commuting with Aq necessarily split in block 
form: Bi = { „„ 1 , therefore the DS in NF will be 



B'/ 



v' =A'ov' + ^ p'.B^v' = (1 + a) AW + ^ . fi^B^v' 

' * (21) 

v" =Kv" + ^ /3fS,r,/' = (1 + a)Kv" + 0^B]v" 

where (3^ and /3" are functions (of A and) of all the constants of motion p{v' , v"). Considering in particular the 
terms /?", the assumption that there are no resonances of the form (20) excludes the occurrence of fractional 
constants of motion in Xyi,, of the form (3" = p'{v')/vk for some k = r + 1, . . . ,n, then v" = solves the 
second set of equations of the system in NF. Transformation into the initial coordinates u shows the final 
statement. 

• 

Remark 2. A quite common situation, which is actually a special case of Lemma 2, is that the n eigenvalues 
(Ti are such that the only analytic or fractional constants of motion p which can be constructed depend only 
on r variables v', with the above notations (e.g. if ai = 2, —1, ^/2, -s/S). In this case, no resonance of the form 
(20) is clearly admitted and Lemma 2 holds true; notice in particular that in this case the NF (21) becomes 
automatically "triangular" (cf. [13]). It is also clear that, given n eigenvalues o-j, it can happen that several 
different reductions are possible: e.g., if tTi — 1, —2, V2, — v^, v^. Lemma 2 can allow a reduction of the NF 
either into a 4-dimensional problem, or equally well into two independent 2-dimensional problems. {} 

Remark 3. We can clearly combine Theorem 1 and Lemma 2 to obtain a resonant bifurcating solution with 
n>2, r = 2, p=l. It can be noted that the leading terms of the resonant bifurcating solution obtained in 
this way can be characterized as the kernel of the linear operator T defined by 

nn d 

acting on the linear space of the vectors w = {wi, . . . ,Wn) where the Wi are formal series in powers of 

Wi= cm.e'"*^"* i = l,...,n (22) 

— — OO 

where, using previously introduced notations (12), 

Sl S2 

It is easily seen, indeed, that the hypothesis on crj is equivalent to the property that T has precisely a 
2— dimensional kernel, generated by wi = e^^^°* = e"^* and W2 = e~^'^^°* = e'^^*; see (14). 

An especially interesting case occurs clearly when one of the eigenvalues is zero. This can be considered 
as a particular case of Lemma 2; we have indeed: 

Corollary 2. (Stationary bifurcation) Assume that for Aq = the matrix Aq of the given DS admits just 
one eigenvalue ( say oi ) equal to zero. Then, with the standard condition 

dAn(A) 



.=0^° (23) 



dA 

there is a stationary bifurcating solution of the form u = ui + . . .. 

Proof. With the notations introduced in the proof of Lemma 2 we put here v' = v\, v" = {v2, . ■ . ,Vn), 
and now v\ G Iaq being fji = 0. Writing the problem in the form (21), there cannot be terms of the form 
P" = p'{vi)/v" (indeed, a2, ■ ■ ■ , (Jn 7^ 0), and therefore v" = solves the second set of equations for v" . The 
remaining equation is then a 1— dimensional equation 

and assumption (23) is easily seen to be equivalent to 9/3*/9A|a=o 0, which ensures the existence of a 
stationary bifurcation, solving vi = with A = \{vi), as in the usual situations. As in the previous cases, 
the proof is completed coming back to the original coordinates u. • 
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5. DS with dimension n > 2: a general result. 



We now consider a n— dimensional DS: according to the above Section, we can assume, for concreteness, 
that there is a resonance involving all the n eigenvalues (see also Remark 4 below). 

Before giving the main result of this Section, the following property may be useful (the proof is straight- 
forward) . 

Lemma 3. Given a DS in NF v = g{v) = Aqv + G{v), the constants of motion p G of the linear 
part are in general not constants of motion of the full DS, but their time dependence can be expressed as a 
function only of the p themselves: {dp/dt)g = ^{p), where {d/dt)g is the Lie derivative along the DS. If the 
DS satisfies Condition A, then these constants of motion p are also constants of motion of the full DS in 
NF V = g{v). 

Theorem 2. Consider the DS (1) and assume that for the value Aq = the eigenvalues ai of Aq are distinct, 
real or purely imaginary, and satisfy a resonance relation (3). Assume also that p = n — 1, i.e. that there 
are n — 1 real parameters X = (Ai, . . . , A„_i), and finally that putting 



(i) _ dA,i{X) 



i = l,...,n: fc = l,...,n-l (24) 

A=o 



the n X n matrix D constructed according to the following definition (notice that only the diagonal terms 
Au{X) of A{X) are involved) is not singular, i.e.: 



det D = det 



(2) 



+ (25) 



Then, there is, in a neighbourhood of uq = 0, Aq = 0, t = 0, a bifurcating solution of the form 

Ui{t) = ( exp(a(A)^oi))^oi(A) + h.o.t. i = 1, . . . , n (26) 
where a{X) is some function of the X 's such that S(A) — > 1 for A — > 0. 

Proof. As in the particular cases examined in the previous Sections, let us consider the given problem 
transformed into NF: v = Aqv + G{v, A) in the new coordinates v. Let us write G{v, A) in the following, 
more convenient form: 

n 

G{v,X)=Y,>^i{\p{v))KiV 

i=l 

where Ki = diag(0, . . . , 1, . . . , 0) is the diagonal matrix with 1 at the i— th position, or also 

Vi = cFiVi + Ki{X, p{v))vi (no sum over i = 1, . . . , n) p{v) e Jaq (27) 

Now recall that the functions p{v) and k{X, p) can be fractional in the components Vi, but in such a way 
that each term KtKiV is a polynomial, so that the only admitted fractional terms have necessarily the form, 

e.g., — — — — —, and so on; the assumption that the ai are distinct ensures that the functions n cannot 

Vi 

be of zero degree in the Vi (i.e. of the form V2/V1, e.g.), then when t; — > all terms KiVi vanish more rapidly 
than V and one finds ^ 

— — I Ki (a, p(?j))tii ) = Ki{X,0) (no sum over i) 
Then, eq. (27) can be written, at the lowest-order 

n-l 

Vi = cFiVi + Ki{X, Q)vi + . . . = aiVi + ^ QikXkVi + ... + ... (28) 

fe=i 
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where qik arc the elements of a constant matrix with n rows and (n — 1) columns. On the other hand, 
considering the original DS 

ii = A{\)u + ... 

its diagonal bilinear terms (in u and A) a^^^XkUi (using definition (24)) are just NF terms and are not changed 

by the normalizing procedure; therefore, a^'^ = q^k. Now according to Condition A, the NF is convergent 
(or better: is obtained by a convergent NT) if one can rewrite (27) in the splitted form as in (10): 

v = Aov + a{\p{v))Aov + Y^*, Pj{X, p{v))BjV (29) 

where a, (3j are suitable combinations of the Hi, and can satisfy the (n — 1) conditions 

(3,{X,p) = (30) 

In fact, the hypothesis (25) ensures precisely that one is able to do this and also to satisfy Pj{X,p) = 0, 
by means of the implicit-function theorem, giving some (n — 1) relations (here the p are considered as 
independent variables) 

X,=X,{p) (31) 

Once these n — 1 conditions are satisfied, i.e. on the manifold defined just by (31) (which is an analytic 
manifold, see [4]), the convergence of the NT taking the initial DS into 

i) = Aqv + a{X)Aov ^ a{X)Aov (32) 

is granted. This DS can be easily solved, giving 

v{t) = ejcp {{a{X)Aot))vo (33) 

with the n — 1 relations 

A, = A, {p{v)) p = p{v{t)) = pivo) e Jao (34) 

The desired result is then obtained, with similar remarks as in the proof of Theorem 1, coming back to the 
initial coordinates by means of the inverse (convergent) transformation v u = v + iIj{v), where ^^{v) are 
series of monomials of the Vi (i = 1, . . . ,n). • 

It is immediately seen that, in particular, condition (13) of Theorem 1 is nothing but a special case of 
(25). As a generalization of Corollary 1, the case of purely imaginary eigenvalues is particularly interesting, 
becaiise it corresponds to the case of coupled oscillators with multiple frequencies and gives, in the above 
hypotheses, the existence of multiple-periodic bifurcating solutions. We have indeed [17]: 

Corollary 3. With the same notations as before, let n = r = 4 and ai = —02 = i, (73 = —oa = mi (with 
m = 2, 3, . . .): then, with X = (Ai, A2, A3) e R^, and det D ^ 0, there is a double-periodic bifurcating solution 
preserving the frequency resonance 

coi : u)2 = i '■ m 

Example 3. Consider a 4-dimensional DS, with u S R^, A € R^, describing two coupled oscillators with 
unperturbed frequencies wi = 1, W2 = 2: 



ii = A{X)u + F{u, A) where A{X) 



/A1 + A3 -(I + A2) Ai \ 

1 + A2 Ai - As Ai 

A2 A3 -2 

V -A2 2 A3/ 



For A = 0, the eigenvalues of Aq are ±i,±2i, and it is easily seen that condition (25) is satisfied (the above 
procedure and notations can be extended without difficulty to the complex space). The NF will have the 
form, denoting by Wa = Vi + iv2, Wb = Vs + ivi the new coordinates, in complex form 

Wa = iWa + (Ai + P^^\X, p) + i{X2 + 4^\A, p)))Wa 

Wb = 2iwb + (A3 + (A, p) + iP^i^ (A, p)) Wb 
10 



where we have put 

A(A,p) = A(A,0)+/3(^'(A,p) 
and the (ii are real functions of the three (functionally independent) constants of motion 

Pi = \wa\'^ =vl + vl = r'l; p2 = Iwftl^ = +d| = ; ps = {wlwb + c.c.) = 2r^r6Cos2<^ 

where f is the time phase-shift between the two components Wa and Wb- Notice that fractional constants of 

motion p{w) may appear in this problem, e.g. iv'^/wi, or w^Wb/wa, etc. (which arc functionally - but not 
polynomially - dependent on the three above), but this would not alter the result, as shown in the proof of 
the theorem. The above NF can be trasformed into the splitted form (29) 

Wa = iWa +ia{X,p{w))wa + Pa{\ p{.w))wa 

Wb = 2iwb + 2ia{X, p{w))wb + l3b{X, p{w))wb 

with 

Pa =Ai + +i{2X2 + 2/3f ) - /3f ^) ; /36 = As + + i(2A2 + 2^'^ - pi'^) 

One can impose the convergence of the NT solving for Xi = Xi{p) the conditions (3a = (3b = 0, which actually 
give three real conditions, and obtain a bifurcating double-periodic solution, with frequencies 

wi = l-|-Q!(A) and a;2 = 2(l -|- Q!(A)) 
Just to give a concrete example, let us imagine that the NF is such that 

Pi — -Pl , P2 — U , Ps — -p2 , P4 — P3 

then the leading terms of the solution, in the original real variables Uj, are 

^^1 = i~a COS ivt ; U2 = Va sin cot ; U3 = rb cos 2co{t + if) ; U4 = rb sin 2u}{t + (p) 
with the constraints 

Ai = ; A3 = ; A2 = r^n cos 2(^ ; w = 1 A2 
producing (see especially the role of A2) a sort of amplitude-phase-frequency locking in the solution. {} 

Rem,ark 4- As already remarked, it can happen that, among the n resonant eigenvalues (7j, as considered 
in Theorem 2, one can find some r < n eigenvalues in such a way that the assumption of Lemma 2 is 
satisfied, and therefore the problem can be reduced - as explained in the above Section - to a r— dimensional 
problem. In this case, if one can also find r — 1 parameters A^ in such a way that the corresponding r x r 
matrix D is not singular, then the existence of another bifurcating solution is ensured by the same Theorem 
2 (see the end of the final example 6 for a - quite simple - case). 

6. Degenerate eigenvalues and the presence of symmetries. 

We now consider the case of multiple eigenvalues of the matrix Aq. This situation is a little bit more 
involved: indeed, the presence in this case of constants of motion p G Iaq of the form p = Vi/vj prevents the 
direct application of the argument used in the previous Theorems; another difficulty is related to the greater 
number of matrices Bi G C{Ao) (if an eigenvalue (jj has multiplicity d, then any d x d matrix acting on the 
subspace of the corresponding eigenvectors clearly commutes with Aq): this would require the presence of a 
greater number of parameters A, in order to satisfy the Condition A. 

However, the presence of degenerate eigenvalues is usually connected to the existence of some symmetry 
property of the problem, and we will restrict to consider this simpler - and probably more realistic and 
physically interesting - case. 
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Wc refer here to the case of "geometric" or Lie point-symmetries [18-19]: a vector function s{u) = 
Lu + S{u) is said to be (the infinitesimal generator of) a symmetry for the given DS if s{u) is not proportional 
to f{u) and the vector fields Xf = / • V and Xg = s - V commute: 

[Xf,Xs]=0 (35) 

or, introducing the Lie-Poisson bracket {•, •} between two vector functions h^^\u), h^^\u) defined by 

{ft(i) , } . = (ftd) . V) ftf - (/i(2) . V) /I W (36) 

if, equivalently, 

{/,s}=0 (35') 

The symmetry is linear if S = 0. Notice that linear and nonlinear symmetries are changed one into each 
other under (nonlinear) coordinate transformations u — > t;. It can also be remarked that a DS in NF always 

admits a linear symmetry [8,11,12]: 

Lemma 4. Any NF admits the linear symmetry sa„ = Aqv: 

{Aov,g{v)} = 

This is in fact a restatement of (7) and (5) using the above definition (36), indeed A{g) = {AqVjQ}. 

We need the following important properties of Lie point- symmetries of a DS. 
Lemma 5. Assume that the given DS admits a symmetry s{u) = Lu + S{u) where L is semisimple and not 
zero. Then one has in particular 

[Ao,L] = 

and there is a NF of the DS which admits the linear symmetry Sl = Lv. 
The proof is well known and can be found, e.g., in [8,11,12]. 

An (unpleasant) consequence of the degeneracy of the eigenvalues CTj of Aq is a larger arbitrarity in the 
choice of the matrices Bi e C{Ao) to express the NF (8): consider e.g. the case in 

= diag(l,l,-2) 

then the same resonant term can be written in two apparently different forms (the notations are obvious): 

x'^yz\ /O 1 

^x^z 

/ yo 

/I 
=xyz 
\0 

Assume now that the symmetry "removes the degeneracy" , that is, recalling that Aq and L commute, 
assume that there is a simultaneous basis of eigenvectors for Aq and L such that any two eigenvectors with 
the same eigenvalue under Aq are distinguished by a different eigenvalue under L. Then, there are precisely 
n independent matrices, which we now denote by Bi, commuting with both Aq and L: 

BiGC{Ao)nC{L) (37) 

These matrices Bi also commute with each other (they in fact can be taken diagonal: this is true upon 
complexification of the space, in the case the eigenvalues of Aq or of L are not real: see Example 4 below). 
We can then express the NF by means of precisely these n matrices S,: 

n 

G{v,X) = Y,l3,{\,p{v))BiV BieC{AQ)r\C{L) 

i=i (38) 

=a(A, p{v))Aqv + ^* (A, p{v))BjV 



p\B\v or also 



P2B2V pi,p2elAo, Bi,B2eC{Ao) 
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using previously introduced notations. The convenience of this choice is immediately evident, indeed, recall- 
ing Lemma 5, the NF admits the symmetry sl = Lv, i.e. 

{Lv,G{v)} = (39) 

which implies 

n 

Y,Lv-VPi BiV = 

i=l 

or - due to the independence of the BiV - 

a{X,piv)), pj{X,p{v)) €Il (j = l,...,n-l) (40) 

where II is the set of the constants of motion of the linear problem v = Lv, and this means that the functions 
Pi in (38) (and the p as well) can be chosen as simultaneous constants of motion of the two linear problems 

V = Aqv and v = Lv: 

p{v), pi{X,p)eiAoniL (38') 

The assumption that L has removed the degeneracy has the other consequence that no constants of motion 
of the form Vi/vj are admitted in Xa„ HIl (and then in the NF (38)), therefore, the functions /3j appearing 
in the ^* in (38) can be written in the form 

/3,(A,p) = /3,(A,0) + /3f^(A,p) (41) 

with fjj^"^ (A, 0) = 0. On the other hand, considering the first-order terms (in u) of the DS in its initial form, 
and writing 

ii =Aou + A^^\X)u + higher order terms 
=Aou + ^ bi{X)BiU + ^ ce{X)Ceu + h.o.t. 



e (42) 
^,bj{X)B,u + Y,' 



^Aqu + a{X)Aou + bj{X)BjU + ce{X)Ceu + h.o.t. 



where ^(^^(O) = 0, the sum includes all linear terms which are not resonant (i.e. [Ao,C^] ^ 0), and 
therefore disappear after the NT. The remaining terms are instead not changed by the NT and therefore one 
gets 

a(A)=a(A,0), 6,(A)=/3,(A,0) (j = 1, . . . ,n - 1) (43) 

Then, assuming that there are n — 1 parameters A^, and that the (n — 1) x (n — 1) matrix D, introduced by 
means of the following definition, is not singular, i.e. that 



~ db 
det£> = det — t" 

oXk 



y^O bj=bj{X)=(3j{X,0) (44) 

A — 



one can proceed exactly as in Theorem 1 and conclude with the existence of a convergent NT and of a 
bifurcating solution on some manifold Xj = Xj{p). 
We can then state: 

Theorem 3. Let u = ^(A)m-I-F(A, u) be a DS with n—1 real parameters Xj; assume that Aq has degenerate 
resonant eigenvalues, and that the DS admits a Lie point- symmetry generated by s{u) = Lu + S{u) where 
L ^ is semisimple, such that L removes the degeneracy of the eigenvalues of Aq . Let Bi be n independent 
matrices in C{Aq) n C{L), and write the linear part (in u) of the DS in the form (4-^). Then, if the matrix D 
defined in (44) is not singular, there is resonant bifurcating solution of the same form (26) as in Theorem 2. 

This result can be well illustrated by two examples (notice that in the second one we will consider also 
the possibility of extending the above procedure in the presence of a discrete symmetry (i.e. not a continuous 
Lie point-symmetry)). 
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Example 4- Introducing the following 2x2 matrices 

J.= {_, ') ,) (45) 

consider the 4— dimensional DS, describing two coupled oscillators with the same unperturbed frequency 
a; = 1, 

■J2 



u = f{u, A) = Aqu + ^ ^i{rl, rl, x n, X)MiU with Aq = ^ " 



J: 



(46) 



where the eigenvalues of Aq are a, = ±i (doubly degenerate), r„ = u\+u\, rj, = u^+ul, Va'^ri, = U1W4 — 
and the Mj are matrices of the form 

J2N" J2N' ) ^" ^'■t'iti'ary 2x2 matrices 

This DS admits the linear symmetry generated by Lu ■ V where 

J2 



J: 



2 



L = 

There are 4 matrices Bj e C{Ao) n C{L), namely 

Ao,L, and/= (^^2 ,iJ= (^^^ 

let us then write down explicitly the linear part (in u) of the DS (46) as in (42), i.e. 
V/|„=o =A(A)=Ao + AW(A) 

4 

=Ao + a{X)Ao + (&i(A)J + 62(A)L + b3{X)H) + ^ q(A)Q 

where a, bj, ce are combinations of the (pi{0, 0, 0, A). When in NF, the last sum disappears, whereas the other 
terms remain unchanged, then the NF will have the form 

v = Aov + aAov+{0^iX,p)I + p2{X,p)L + P3{X,p)H)v (47) 

where p{v) = vf + + vj + v^, V1V3 + V2V4, G Iaq H II, and 

(3j{X,0)=bj{X) i = 1,2,3 

The assumption (44) dct{dbj /dXk)\\=o 7^ then ensures, as already discussed, the existence of a resonant 
bifurcating solution, which in this case is a periodic solution with frequency w = 1 + . . . and preserving the 

strict frequency resonance 1:1. <J) 

Example 5. As a even simpler example, let us consider the case of a 3— dimensional DS with real eigenvalues: 
let 

ii = Aqu + 'Pi{r^, z, X)Lu + (p2{r'^,z, X)Iu (48) 
where u G R^, X G R"^ , r'^ = ul + U2, z = U3 and, using the 2x2 matrices introduced in (45), 



-2 7 V 
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This DS admits the Hncar symmetry generated by Lu ■ V. There arc 3 matrices Bi in C{Ao) n C(L), 
namely Ao,L,I (notice that, actuaUy, in this case, C{Ao) C C{L), but the DS (48) is not in NF, because 
and z are not in Iaq ) • Proceeding as before, the NF will have the form 

it = Aqv + aAov + /3iLv + (3-2lv (49) 

where now j3j (and a, of course) are functions oi p = {v1 + (y'D^y's and A only. Then we need two real 
parameters Afe, and - with the assumption (44) - a bifurcating solution is obtained 

Ml =Mio exp(l + Q)t + h.o.t. 
U2 =U2o exp(l + a)t + li.o.t. 
U3 =W3oexp ( — 2(1 + a))t + h.o.t. 

along with some conditions 

Ai=Ai(r2^); X2 = X2{r^z) 



This example can be useful to show that also the presence of discrete symmetries (e.g. exchange or 
reflection symmetries) can be of some help in this approach: a DS admits a discrete symmetry R , where R 
is a nonsingular matrix, if 

fiRu) = Rf{u) (50) 

and it is known that the NF also admits the same symmetry [5] . The possible role of this fact in our argument 
can be illustrated by the following modification of the above example. 

Example 5 '. The same as example 5, but now assume the DS admits the discrete symmetry 

Then, in both the initial DS (48) and its NF (49), the term containing the matrix L disappears; therefore. 
Condition A requires the vanishing of only one term in the NF and a bifurcating solution can be obtained 
with the presence of just one real parameter A. 



7. Final remark on the role of Condition A. 

In order to stress and illustrate the relevance of the role played in our argument by Condition A, let 
us consider this final example, which is essentially an adjustment (in view of the present discussion) of an 
example given, with quite different purposes, in [20]. 
Example 6. With u e and A e R^, let 

/-(r2 + 3z2)Mi\ /A1+A2 1 \ 

u = A{X)u+ [ -{P + 3z^}u2 \ where A{X) = I -1 Ai (51) 
V {3r^+z^)z J \ -XiJ 

with = Ml + z = M3. It is not difficult to see (cf. [20]) that, if A2 = 0, this system admits a family of 
heteroclinic orbits connecting the origin to the circle = Ai and another family of heteroclinic orbits, living 
on the manifold r'^ + z'^ = Ai, connecting the circle = Ai to the point P = (0,0, Ai). When A2 7^ 0, this 
heteroclinic structure breaks down, and an application of the Melnikov theory [20-23] shows the occurrence 
of transversal intersections of stable and unstable manifolds, with the consequent appearance of the chaotic 
behaviour described by the classical Birkhoff-Smale horseshoe-like structure (actually, we need here a simple 
3-dimensional version of the standard Melnikov theorem: see e.g. [20]). 
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On the other hand, the NF of the above DS (51) exhibits a perfcietly regular (i.e., non-ehaotic) behaviour: 
indeed, the NF, according to Lemma 4, must possess the hnear symmetry generated by Aq, which imphes 
that the NF is symmetric under rotations around the axis; then, the NF is essentially a 2-dimensional 
problem and therefore no chaos is admitted. This clearly implies that the NT cannot be convergent. 

If one now imposes Condition A on the NF of the DS (51), it can be easily seen that this condition is 
satisfied only along the circle defined by 

Ai = 3{vf + vj) +v'l Aa = 4{v3 - vf - vj) 

where in fact a (Hopf-type) periodic solution occurs. Then, in conclusion, we are here in the presence of a 
regular solution (where Condition A is indeed satisfied), which is completely surrounded by chaotic solutions: 
this clearly confirms the crucial role played by Condition A in the argument. 

Just for completeness, and in agreement with Lemma 2 and Corollary 2, let us remark that this example 
admits (quite trivially) a reduction according to Remark 4: indeed, in correspondence to the eigenvalue 
(73 = 0, we also get the stationary bifurcating solution Ai = z"^ with ui = U2 = 0. {} 
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